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f^ l' Abstract 

'^'h ' We use the master equation approach to calculate the energy level shifts of an atom in the 

|/-\ I presence of a general dielectric semi- infinite medium characterized by a dielectric constant e(tj). 
Particularly, we analyze the case of a non-dispersive medium for which we obtain a general 

T-H ' expression for the interaction as well as the asymptotic behaviors for koz <^ 1 (non-retarded 

^ , regime) and koz S> 1 (retarded regime), where loq = fcoc is the main transition frequency of the 

^T ' atom. The limiting cases e ~ 1 and e S> 1 are discussed for both retarded and non-retarded 

C^ ' limits. For the retarded limit, we compute the non-additivity contribution of van der Waals 

^—r^ , forces. 
"^ . 

o. 

^ ! 1 Introduction 

i-Q I Experiments performed by Verwey and Overbeek [^ with colloidal suspensions during the years 

♦"•■ 1940-45 suggested that London's |21 interaction between two polarizable atoms (that falls as 1/r^) 

■y I was not correct for large distances. In order to get an agreement between theory and experimental 

t^ ■ data the van der Waals force should fall faster than 1/r^ for large distances. The authors conjectured 

^ I that this change in the behavior of the force was due to the retardation effects of the electromagnetic 

(sJ '■ interaction or, in other words, due to the finiteness of the velocity of light. 

^ I Motivated by the discrepancies pointed out by Verwey and Overbeek, in 1948 Casimir and 

Polder 15 considered for the first time the influence of retardation effects on the van der Waals 
forces between two atoms and between an atom and a perfectly conducting wall. The non-retarded 
j^ ■ dispersive van der Walls force between a polarizable atom and a perfectly conducting wall can be 

explained with the aid of the image method ^, which shows that the interaction potential varies 
as 1/z^, where z is the distance between the atom and the wall. When retardation effects are taken 
into account, Casimir and Polder [3j showed using perturbative QED that the interaction potential 
becomes proportional to 1/z*. Since Casimir and Polder's paper, forces between atoms or molecules 
and any kind of walls are usually called Casimir- Polder forces. 

In 1956, Lifshitz and collaborators [H] developed a general theory of van der Waals forces. They 
derived a powerful expression for the force at finite temperature between two semi-infinite disper- 
sive media characterized by well defined dielectric constants and separated by a slab of any other 
dispersive medium. Several results have been predicted with the aid of this formula, as for example, 
the variation of the thickness of thin superfluid helium films in a remarkable agreement with the 
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experiments -S] . The Casimir-Polder force can also be obtained as a limiting case of Lifshitz formula 
when one of the media is sufficiently dilute such that the force between the slabs may be obtained 
by direct integration of single atom- wall interactions [3 . 

Since then, a lot of work has been done on van der Waals and Casimir-Polder interactions. Good 
reviews have been written on dispersive van der Waals interaction [S| |51 110| and many elaborated 
papers concerning level-shifts near surfaces have appeared, as for instance, J11II12I 1011141 . to mention 
just a few. It is worth mentioning that Casimir-Polder forces have been observed experimentally 
[TS] . Recently, the influence of real conditions on the Casimir-Polder interaction has been considered 
|16j . Further, higher multipole corrections ^, roughness and corrugation of surfaces [181 ri]?| . 
the influence of the Casimir-Polder interaction on Bose- Einstein condensates [301 12] and possible 
applications to nanotubes |22j are some of the many branches of great activity on this subject 
nowadays. 

Here we shall compute the van der Waals interaction between an atom and a dielectric semi- 
infinite medium characterized by a well defined dielectric constant e{uj). We shall employ the master 
equation approach 23 recently applied by the authors |211 to this kind of calculation. For the 
particular case of a non-dispersive medium (e (w) — e — cte for any frequency), we will also compute 
in a closed form the variation of the interaction potential with e in London-van der Waals and 
Casimir-Polder limits. A measure of the non-additivity of van der Waals forces will be given in the 
large distance regime (Casimir-Polder limit). 

2 General Level Shifts 



In a previous paper [21] we analysed the interaction of an atom with a perfectly conducting wall 
starting from the general expressions for the energy level shifts of a small system interacting with 
a large one considered as a reservoir. These expressions have been derived from so called master 
equation |25| . Let us briefly sumarize the formalism to be used. 

Consider a system S + TZ, where by assumption 5 is a small system and TZ a reservoir. Starting 
from an interacting Hamiltonian of the form 



vit) = -J2sAmAt) 



(1) 



, where Sj(t) and Rj(t) are S and TZ observables, the master equation comes from the perturbative 
expansion of the density matrix of the total system S -\- TZ. Working with a second order pertur- 
bation theory, one can take the reduced trace over the Hilbert space of TZ and do the Markovian 
approximation. This approximation is due to the existence of two very different time scales, namely: 
the correlation time of the fluctuations of the TZ variables and the characteristic time evolution of 5, 
which is too large compared to the former. The final equation describes the evolution of the density 
matrix of S, made up of a free term plus a linear term that accounts for the coupling with 7?, |25| . 
The master equation is given by: 
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where gjk (r) = Trj^; [prRj (t) R^ (0)], Pr[s) — ^^S[R) {Ps ® Pr) is the reduced density matrix of 
TZ{S), pfj — {a\ps\b), S^f^ — {a\Sj\b), hcUab = Ea — Ef, and \a) represents an energy eigenstate of 
S with eigenvalue Ea- In order to obtain the level shifts of S it is necessary to consider only the 
non-diagonal terms which are coupled with themselves or, in other words, the terms in which the 
indexes c = a and d — b. For simplicity, one assume that levels \a) and \b) arc non-degenerate. 
Hence, from eq.Q, one may write: 

Jahab = -Fafc - iAab ; ^a ^ T }_^ P^. }^ }_^V _ _ (4) 

"- p „ n ^t^^ ^"^ ^^ ^" 

where A^fc = A^ — A^, Tab is a damping, V is the Cauchy principal value, p^ is the statistical weight 
of the state \p) which is an eigenstate of TZ with energy eigenvalue E^ and V is defined by eq.l^. 
Using the definitions of the linear susceptibility and symmetric correlation function for both small 
system and reservoir, the level shifts can be splited into two terms: the first one is due to polarization 
of the system induced by the fluctuations of reservoir (/r) and the other is the contribution of the 
reservoir reaction [rr) over the system due to the fluctuations of the system induced by the reservoir. 
The expressions are: 
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where &Ea = hAa = ^A^J + ^AY is the energy level shift of S in the energy eigenstate \a) with 
eigenvalue Ea, Xjfc (^) is the real part of the susceptibility of TZ and C^ (w) is its symmetric corre- 
lation function, C,^"' {uo) and Xjfe ° (^) are the equivalents of Cj| (w) and x'll ('^) ^^^ ^ i'^ ^^^ state 
I a) and -R^^ = {ij\Rj\v), where \p) is the energy eigenstate of TZ with an eigenvalue E^. 

Let us consider the electromagnectic field as the reservoir and the atom as the small system. 
This is allways possible when the correlation in time of the field observables is very sharp compared 
to the correlation of the atom observables: this last correlation width is ^ 1/loq while the field 
correlation is '^ ^/^jJc, where loq is the main frequency transition of the atom and lJc is an arbitrarily 
large cut-off frequency. 

In the dipole approximation the interacting Hamiltonian may be written as: 

T/(x,t) = -d(t).E(x,i) (11) 

where d (i) is the dipole moment (the system observable) of the atom induced by the electric field 
(the reservoir obervable). Since this Hamiltonian is bilinear in the atom and field operators, one 
may use expressions ^ and ((HJl to calculate the interaction. 



Consider a field mode as 



EkA (x, t) ^ FkA (x, t) al^ + h. 



(12) 

(this is always possible in a linear medium) , where FkA (x, t) is a function of time and position 
that takes into account all sourced contributions and boundary conditions imposed to the field, and 
a\^y^ is the creation operator of a photon with wave-vector k and polarization A which satisfies the 



commutation relations akAi^k/A' ~ '5kk"^AA'- Consider also the atom as an isotropic two-level 
system with transition frequency wq — fcoc = (i?e — Eg)/h where E'e and Eg are the energy of its 
unperturbed excited |e) and fundamental \g) states respectively. This implies that all non-diagonal 
terms of its susceptibility and second order symmetric correlation function given by © and Q), 
respectively, are identically zero or, in other w! ords: Xjk ^^) ~ X" {^) ^jk and C-^"" {to) = C"" {uj) Sjk- 
Then, if the atom is in its ground-state, equations H7I1U|) and H12|l will lead to: 
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where uq = 2\deg\ /3?iujo is the static polarizability, k = |k| and (nkA) is the statistical average 
number of photons in a given mode. Combining equations (|13I16(I with jSJ and lO, one can write 
the level shift SE = HSluo = SE^'' + SE""" as: 
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In equation (|18|l 0;+ (fc) is the well known approximation for the polarizability of a two-level 
system with a negligible natural line- width r„, namely cjq S> r„. It is formed by a resonant part, 
V [(fc — fco)~^] , and a non-resonant one, V [(fc -f fco)^^] ■ Then, the dipole moment d^ of system S 
induced by the field is given by d^ (x, t) — X^kA '^+ (^) -^kA (x, t). 

By the same argument, a_ (fc) may be interpreted as a field polarizability related to the reaction 
on itself of its influence on the system S. Equivalently, one may talk about the "dipole moment" 
A-ji of the field, that can be written as A-ji (x, t) — X^kA ^~ i^) -^kA (x, t). 

For the vacuum state of the field we have (rikx) — and the total energy level-shift may be 
written as: 

2 ^ k + k, 

where only the non-resonant term contributes. In the following section we will apply the previ- 
ous equation to find a general expression for the interaction between an atom and a semi-infinity 
dispersive dielectric medium for all distance regimes and all possible values of the dielectric constant. 



3 Force between the atom and the dielectric medium 

Let us consider a semi-infinite dielectric medium with a dielectric constant e (fc) defined in the region 
z < 0. Let us also consider that there exists an atom at a position z > 0. For this simple system, it 
is possible to separate explicitly the modes of the electromagnetic field f5H|. Hence, for the electric 
field one may write |S]: 



E 



/2nhck\ 



1/2 



kA 



(.,(x,i) = -M-^j AkA,<,,(x)e^'^^*ai,,^^^^(0) + /i.c, 



AkAi (x) 
AkA2 (x) 

fc3 
k 



1 

1 

kV2 



e ^■''ekA 



e ^■''ekA X 



fc.C. - kn 



3 1 2ik3Z 



fcg 

e (fc) fca 



fc2 

k^ + fcgZ 






k' — kj^ — ^32 and k^ — kix + k2y 



U2 1,2 . 

- — «,_L , 



(21) 
(22) 
(23) 

(24) 
(25) 



As a consequence, the sum in polarizations of the modulus square of the amplitude of the field is 
given by: 
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Using the previous equation on (|20|l . remembering that ^j. -^ V/ {2tt) Jk^dkjdfl in spherical 
coordinates and keeping only the z-dependent part of 5E^ which is the only one that contributes to 
the interaction of the atom with the wall, one obtains after some mathematical manipulations: 
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where the principle part V was omited because there is no poles, since k — |k| and ko are allways 
positive and the integral in k is over on the positive real axis. 

Equations (|27I28() can be considered one of the main results of the present paper. They give the 
interaction between an atom and a semi-infinite dielectric medium for all distance regimes and any 
frequency dependence of the dielectric constant. As an important particular case we shall consider 
in a moment a non-dispersive medium. It is worth mentioning that Nha and Jhe ^14, have also 
considered non-dispersive media. They calculated the level shifts of an atom between two parallel 
dielectric surfaces using linear response theory |27[l28j . 

Similar results involving dispersive media have been deduced in the literature directly from the 
interatomic potentials or from Lifshitz formula (see [H] and references in it). However, as far we 
know, in all of them it is necessary to consider at least one rarefied medium and, if necessary, add 
perturbatively the corrections due to the non-additivity of the van der Waals forces, which arises 



from many-body interatomic potentials [2211201 (for a simple discussion see |2I])- Note that our result 
gives a very compact formula for the interaction between the atom and the semi-infinite dielectric 
medium which includes all many-body interactions. 

In order to proceed, it is necessary to have the explicity form of e (k). Let us analyze the particular 
case of a non-dispersive medium, where the dielectric constant may be considered independent of 
/c in a large range of frequencies or, in other words: e (fc) ~ e (0) = e. For this case, one can make 
firstly the integration on k in equation (|27|l . which leads to 
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where 7 is the Euler-Mascheroni constant and xq — 2fcoz. Integrating equation (|29|) by parts, one 
may write: 
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where /', /", ..., mean first, second, ..., partial derivatives of / respect to t. In order to check the 
self-consistency of our results, let us consider the particular cases e = I and e ^ 00. By direct 
inspection of equation lpH|l , we get 
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^From equations (|36|l and H29|l we see that for e = 1, the interaction potential is identically zero, as 
expected, since e = 1 means that there is no medium at all. 

For the case e — > 00, the last term on the r.h.s. of equation 135|l vanishes, since from l|37|l 
/'" (t, e) = 0. The remaining terms give a finite value so that the interaction potential may be 
written as: 
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in perfect agreement with a previous result |24| . Using the hmiting expressions oi J-' (x) and Q {x) = 
T' (x) for smaU and large values of x: 
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one may write for the asymptotic behaviors of Vb {z): 
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which are the expected results for a perfectly conducting wall. 

For any finite value of e it is not possible to simplify expression H35|l. as we did for the e ^ oo case 
(see equation (j38|l ) . However, one can obtain the e-dependence in a closed form for the asymptotic 
expressions for large {xq ^ 1) and small {xq <C 1) distances. In these cases, it is possible to make 
the factorization: V {z, e) — 5 (e) Vb(z), where g{e) is a function of e that will be different for each 
distance limit. Directly from (|35|l . we obtain 
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In order to evaluate the above integrals, it is convenient to define k := e — 1 and then consider 
separately two cases, namely, k <^ 1 and n very large. In the former, we expand f"'{\t\,K) in a 
power series of k, while in the latter we expand /"'(|i|, k) in a power series of 1/k. These series are 
given, respectively, by: 
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The above expressions are strongly divergent at i = 0, which belongs to the interval of integration in 
both integrals appearing in equations (|43|l and H44() . In order to circumvent this problem and extract 
a finite result from these ill defined integrals, we replace t hy t —^ t ± iS and, after performing the 
integrations, we take (5^0. Adopting this regularization prescription, we may write directly from 
equation H43|) an expression for the interaction between the atom and the semi-infinite dielectric 
medium for the case xq <C 1: 
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We should emphasize that the previous equation is valid for e > 1. In order to give a better 
understanding of this result, let us consider a sufficiently dilute medium so that the non-additivity of 
the van der Waals forces may be ignored. In this case, as a first approximation, one can use London's 



result for the interactfon between two identical electrically polarizble atoms, Vab = —^hcuoaQ / Ar^ , 
to make a pairwise integration between the atom and the all other atoms that constitute the semi- 
infinite dielectric medium. Using the Clausius-Mosotti relation 32 , iirNaQ = 3(e— l)/(e-|-2), where 
N is the number of atoms per unit volume of the medium, and assuming for a dilute medium that 
e « 1, one obtains exactly half the value given by l|45() in first order in e — 1. This discrepancy is 
due to the fact that the London's result is valid only for distances smaller! than the wavelength of 
the dominant atomic transition and then can not be used for integration over all possible distances. 
For the opposite limit, xq 3> 1, the expression for g{e) is not so simple like that written in H45|l . 
However, it is not difficult to obtain series expansions in 1/k or in k. For large values of k, we may 
write from H44() : 

which goes to result H42|l when e — > oo. For small values of k, it is possible to write: 
,,, , 23 aohc ( 169k 2263^2 \ 

Let us look at the previous equation more carefully. Following the same procedure as that adopted 
to obtain equation H45() , we can compute the interaction from the direct pairwise integration of the 
atom-atom potential in the Casimir-Polder limit: Vab (r) = —2ifica^/ATTr^. Setting e ~ 1 and using 
again the Clausius-Mosotti relation, we may write for the integrated potential \^(™*' (z, k): 

One can see from above that the leading terms of equations H47|) and (|48|l are the same. This 
means that in first order on k the interaction is due only to direct integration of the atom-atom 
potential in the Casimir-Polder limit. The higher order terms, however, are different because of the 
non-additivity of van der Walls forces so that the integration of a two-body potential can not give 
the exact result. A measure of this discrepancy may be given by the quantity: 

F(z,K)-y(™*)(z,K) 185 303113 o 1325388223 „ 

— ——^ — ^ — LJ_^ ^ kH K^ K^ + ... (49) 

V{z,k) 966 3732624 39662862624 ^ ' 

which essentially provides an average contribution of the many-body potentials to the interaction; 
this contribution is second order on n. One may note also that differently of two-body potential 
that is always attractive, in this case the many-body contribution has a repulsive character (the 
many-body interaction may be repulsive depending on the geometry considered \S3\). 

4 Conclusions 

In this work we studied the interacton of an atom with a semi-infinite dielectric medium characterized 
by a dielectric constant e(w) which, in principle, can depend on the frequency. We used the general 
expressions for the energy level shifts of a small system coupled to a reservoir given by the master 
equation approach. For simplicity, we considered the atom (interpreted as the small system) as a two- 
level system and the radiation field (interpreted as the reservoir) in its vacuum state but submitted 
to the boundary conditions imposed by the presence of the semi-infinite dielectric medium. With 
these assumptions, we were able to establish a quite general expression for the interaction of the 



atom with the semi- infinite dielectric medium given by equations H27I28|I . It is worth emphasizing 

that this resuh takes into account dispersion and is valid for all distance regimes. 

We also analyzed the particular cases of a non-dispersive medium, where we obtained simple 
formulas for the interaction potential, specially in the retarded and non-retarded regimes. For these 
cases, it was possible to write the potential as a product of a function of e only and the potential for 
a perfectly conducting medium in the corresponding limits. For the non-retarded limit, we obtained 
a simple analytical expression for the atom-medium interaction valid for any e. For the retarded 
limit, however, simple expressions were obtained only for small values of k = e — 1 and large values 
of K. Note that the first term in expansion (46), valid for large k and large distances coincides, as 
expected, with the Casimir and Polder result 3 . Further, for small k and retarded regime, we also 
computed the contribution of the many-body potentials to the interaction providing a measure of 
the non-additivity of the van de! r Waals interactions. 
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